New mechanisms of droplet coarsening in phase-separating fluid mixtures 
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We propose here a new mechanism of droplet coarsening in phase-separating fluid mixtures. In 
contrast to the conventional understanding that there are no interactions between droplets in the 
late stage of spinodal decomposition, we demonstrate the existence of interactions between droplets 
that is caused by the coupling between diffuse concentration change around droplets. We show 
the possibility that this mechanism plays an important role in droplet phase separation together 
with Brownian-coagulation mechanism. We also discuss the coupling between hydrodynamic and 
diffusion modes, namely, "collision-induced collision" phenomena. 



PACS numbers:64.70.-p, 64.75.+g, 05.70.Fh 

Droplet patterns are generally observed in the late 
stage of phase separation in binary mixtures of solids 
and liquids at off-symmetric compositions |QJ2). There 
are two well-known coarsening mechanisms for droplet 
phase separation: (i) evaporation-condensation mecha- 
nism Q,^] and (ii) Brownian-coagulation mechanism due 
to collision between droplets that is caused by the free, 
thermal diffusion of droplets without any interaction 
jl],(i]j|]. The former is known as Lifshitz-Slyozov- Wagner 
(LSW) mechanism and the latter as Binder-Stauffer (BS) 
mechanism. Both mechanisms give the same form of the 
power law of for the evolution of the mean droplet radius 
R as a function of the phase-separation time t, in three 
dimensions [fd^fl0|: 



R 3 = k d (k B T/5irri)t, 



(1) 



where kd is a constant, T is the temperature, r\ is the vis- 
cosity, and is Boltzmann's constant. For LSW mech- 
anism kd = 0.053, while for BS mechanism kd = 6$^ 
(&d'- the volume fraction of the minority phase) 
The droplet coarsening behavior in phase-separating fluid 
mixtures has mostly been explained in terms of these 
LSW and BS mechanisms so far. Recently, however, 
some reserachers start to think that this well-accepted 
view needs some essential modifications in physics on the 
basis of experimental results and theoretical considera- 
tions. There are two important physical factors that have 
so far been ignored: (i) interdroplet interactions 1 1 1 - 1E|| 
and (ii) hydrodynamic interactions 10,0 ■ 

In this Communication, we consider the interdroplet 
interactions between droplets and the resulting hydrody- 
namic motion of droplets |L7|]. We think that the con- 
ventional picture that a droplet moves around freely by 
Brownian motion and accidentally collides with another 
droplet is not correct in the exact sense, at least, for a 
system of a high droplet density p^Jl^ |. 

The kinetic equations for classical binary fluids are |l]] 



dcj) 
Of 



= -V 



L V 



,6((3H) 



dv - 



(3) 



where 4> is the composition, v is the velocity, p is the 
density, p is a part of the pressure that is determined to 
satisfy the incompressibility condition, 9 and Tare ther- 
mal noises, and (3 = \/k B T. Here H is the Ginzburg- 
Landau-typc Hamiltonian: 



H = k B T J dr[f(cj)) 



K 



(V0) 2 



(4) 



with f(<f>) = -§</> 2 + f </> 4 . In Eq. (§), F is the ther- 
modynamic force density acting on the fluid due to the 
fluctuations of the composition <f> that is given by 

= -cfNp = -V[?r - KV{cj)V 2 (j))} - KV 2 ct)V(t), 



knT 



where fj, = S((3H)/d<p is the chemical potential and 
7r is the osmotic pressure given by ir = <pdf/d(fi — f. 
The conserved part of the above thermodynamic force, 
F^/k B T — — V[7r — KV (4N 2 (j>)], cannot produce any ve- 
locity fields under the incompressibility condition V • v = 
and it has to be completely balanced with a part of 
Vp. Accordingly, only a part of the other force, FY 1 * = 
~KV 2 (t>V4>, can produce the velocity fields through Eq. 
(^). Using the steady-state approximation (dv/dt ~ 0), 
thus, we obtain the following expression of v p|pq|: 



(5) 



v= / dr'T{f-r')-Fl nt (r'), 



+ i 



(2) 



where T is the so-called Oseen tensor. 

This expression gives us the physical insight on hydro- 
dynamic effects on domain coarsening as follows |^,|l9|,^0| . 
After the formation of a sharp interface, it has so far 
been assumed that the interface profile can be described 
by <pi n t = 4>e tanh(C/v / 2^) where <j> e (= y/r/u) is the 
equilibrium composition, £ is the correlation length. £ is 
the distance from the interface defined by £ = rt-(r—fi n t), 
where r int is a point on the interface and n is the unit 



f 



normal vector at the point ?i n t toward the domain with 
a positive value of <j>. Then, the thermodynamic force 
density FY 1 * can be expressed by 



(6) 



Note that V • ft is the curvature at ri nt and V • ft = 
1/Ri + I/R2, where l/Ri and I/R2 are the two princi- 
pal curvatures of the interface. Here we use the relations: 



V0 = (d<f> int /d()n, 



d 2 4> 



int , / (-'(pint > ^ _, 



(7) 
(8) 



Thus, we obtain the following equation by putting Eq. 
© into Eq. (|) and using k B TK(d(f> mt /d0 2 = oS(Q: 

-Vp'-a(-^ + -^(C)n + ?7V 2 « = 0, (9) 

where a is the interface tension. It should be noted that 
Eq. (^) is dependent on $d through the characteristic 
curvatures of the interface. Under the incompressibility 
condition, thus, the domain geometry (ft pattern) domi- 
nates the coarsening mechanism: For a nearly symmetric 
composition (§d ~ 1/2), we have a bicontinuous pattern; 
and, thus, the second term in Eq. (^) produces the ve- 
locity fields of v ~ ajr\ that lead to the growth law of 
R ~ (a/n)t (Siggia's mechanism). For an off-symmetric 
composition ($d 7^ 1/2), on the other hand, we have a 
droplet pattern; and, thus, the second term in Eq. (|]) is 
balanced with Vp' to satisfy V • v = 0. Thus, there is a 
pressure difference of 2a /R across the interface (Laplace's 
law), while there are no macroscopic velocity fields v = 
except during the coalescence process. The latter fact 
is the basis of Brownian-coagulation mechanism, which 
assumes that there is no interparticle interaction and the 
droplet motion is purely driven by thermal velocity noises 
pin Eq. (|)]. 

Here we would like to raise a fundamental question 
of whether droplets can feel neighboring droplets or 
not. The above argument tells us that the answer is 
"no". However, we point out that in the above argu- 
ment there are deficiencies that can be important partic- 
ularly in phase separation of off-symmetric binary mix- 
ture; namely, (i) the gradual concentration change around 
a droplet is neglected above and the concentration profile 
is wrongly assumed to be always described by 4>int along 
the direction perpendicular to the interface and (ii) hy- 
drodynamic flow induced by droplet collision is neglected. 
We focus here mainly on the problem (i). Since the 
matrix can never be in equilibrium with droplets with 
different sizes at the same time, a droplet usually has 
long-range concentration (diffusion) fields around it, as 
recently pointed out by us [|l2] |l4| and Karpov We 
think that this fact plays an important role in droplet 



coarsening, especially for fluid mixtures: The finite cur- 
vature of droplets naturally leads to the diffuse concen- 
tration change around them and the coupling between 
these concentration fields causes interactions between 
droplets. Experimental evidence of such interactions be- 
tween droplets has recently been reported for fluid mix- 
tures by Aver'yanov and Kalmykov |ll| and Tanaka Jl2| . 

Before discussing the dynamics, we need to know the 
quasi-stationary concentration profile around a droplet 
for estimating the strength of the interdroplet inter- 
action. After the formation of a sharp interface, for 
both spinodal decomposition (SD) and nucleation and 
growth (NG), we have dS<p/dt = L V 2 ^ neglect- 
ing collision effects. Under the quasi-stationary approx- 
imation, thus, the concentration profile around a three- 
dimensional spherical droplet can be obtained by solving 



the Laplace ec 



equation, |,|,|f|,||,|| 



uation V 2 o 



under the Gibbs-Thomson 
= 4>+(4>b — 4>)R/r, where <p is 
the far-field concentration and 4>b is the concentration at 
the droplet boundary, fa is given by the following Gibbs- 
Thomson equation: (j) b = —<f) e + (a/2kQT(f> e T)(l/R). 
Note that the this 1/r dependence of 4>(r), or /x(r), orig- 
inates from the conserved nature of the order parameter. 
For 8<fi (= (j> — <j>), thus, we have the following relations: 



(r) = 5(f) b - (r > R), 5<f>{r) 
r 



b b (r < R), (10) 



where 8cf>b = 4>b — 4>- 

This diffuse concentration change around droplets, 
8cj>(r), qualitatively alters the preceding conclusion that 
there are no velocity fields for droplet phase separation 
except for thermal noises. In the preceding discussion, 
the relation V0 = (dfant / d()ft are used. This relation 
is reasonable as the first approximation and can be used 
for the case of bicontinuous phase separation, where the 
energy associated with a sharp interface governs the fluid 
dynamics. For the case of droplet phase separation, F^ nt 
given by Eq. (^) does not produce any velocity fields, 
as described before; and, thus, the first-order corrections 
do play important roles. Accordingly, we should use the 
following relation that include the first-order term of Sip, 
instead of Eq. (R): 



V0-(- 



'Hnt 



)n + V5<j>. 



(11) 



Here we assume that the phase having a positive value 
of <j> is the minority phase. In such a case, we have the 
relation Ri = R 2 = —R for a spherical droplet. Equa- 
tion (|9|) should be modified using Eq. ( p"l| ) to include the 
additional thermodynamic force 



F, 



grad 



-K\ 



K 



6<k]6(OV6<l>. 



(12) 



2 



Here we use d(j>i nt /d( = 2</> e <5(£). We neglect the second 
term since it is smaller than the first term by a factor of 
~ £/R. This interface force produces the following new 
velocity fields for droplet phase separation in addition to 
the thermal velocity noises T: 



J dr'T{f 



F n , 



(r'). 



(13) 



This is the theoretical basis of fluid motion induced by 
diffuse concentration change around droplets. 

Next we estimate the interdroplet interaction via 
Fgrad- First we consider the coupling between a single 
spherical particle with a radius of R at the origin and 
the concentration field of a constant gradient given by 
E ■ r. This gradient induces the anisotropic components 
in the force given by Eq. (|T2|). Using the analogy be- 
tween this problem and the change of electric potential 
caused by a metallic sphere in an electric field [ fi"5|]2^ ] , we 
can estimate the resulting matrix concentration profile 
around the spherical droplet, 84>e, a s 



;(r) = 54>{r) +E-r-E- rR 3 /r 



(14) 



Using Eq. (|l4|), we obtain the total force acting on the 
droplet by integrating the terms containing E in Eq. (jlT 

as 



F 



drop 



4K(f> e 



k B T R 

= WirKR^E 



dr5(\r\-R)V8(f)E(r) 



Using this relation, we can straightforwardly calculate 
the total force on a droplet at r[ due to the concentration 
fields of a neighboring droplet at by simply replacing 
the concentration gradient of E by V<50(r) (r = rj — r5) 
as 



-16irk B TKR 2 



(16) 



Here we assume the radii of both droplets are R. This 
force can be easily estimated for the case of <f> = — 4> e and 
it gives the strongest force [||: F^ ax = -~/k B TR^, 
where 7 = 16irao — 5 ~ 10. Here, we use 8<pb = 
(2/3)0 e (e/R), e = K/2t, and a = 4K^> 2 J3^ ~ 
cofeB? 1 /^ 2 (co = 0.1 ~ 0.2 is the universal constant). 
Here it should be noted that this force is not due to 
the interaction between two droplets in the exact sense, 
but due to the coupling between a droplet interface and 
the concentration gradient created by the neighboring 
droplet. 

For a fluid system, the resulting droplet velocity v^ ax 
can be estimated by balancing the frictional force F$ — 
-5ttt]Rv r and F% ax , as 



rpmax 
-vnax _ r R 
V R 



^DrRt 



5iri]R 



knT 



(17) 



Here Dr is the diffusion constant of a droplet of radius 
R and Dr = k B T / '(5%r)R) . It is worth noting that a 
simple dimensional analysis of Eq. ( JTt| ) tells us that this 
new mechanism alone leads to the domain growth law of 
R~t x ' 3 . 

On the basis of the above results, we here estimate 
the maximum coarsening rate including both Brownian- 
coagulation mechanism (thermal noise effects) and the 
above diffusion-coupling mechanism (see, on the calcula- 
tion technique, Refs. p|, |l5| , p5[ ) . Let us consider a droplet 
fixed at the origin in the uniform cloud containing n 
droplets per volume. We suppose the coalescence process 
is much quicker than the collision interval, which is a rea- 
sonable approximation for a fluid system Q . The bound- 
ary conditions for n are (i) n — at r = R, t > 0, and (ii) 
71 = n as r — > 00. The flow of droplets to the origin I is 
given by 7(r) = 4:Trr 2 j(r) with j(r) = —DRdn/dr + uvr. 
Note that I(r) should be constant with r. Using I{f) = I 
and the above boundary conditions, we finally get 



— = 2nl = 
dt 



8nn 2 D R Rj 
- exp(- 7 /2) : 



(18) 



where the factor 2 accounts for the fact that the droplet 
at the origin is free to move. Using the relation n — 
$ d /(lfi? 3 ), we obtain from Eq. @ the following do- 
main growth law: 



R=[6^ d (D R R)] 1 /h 1 ^. 



(19) 



(15) 

The corresponding scaled coarsening law is given by 



(R/0 3 = Gl^dT, 



(20) 



where r = t/r ? [t 6 = i 2 /D^ (D 6 = k B T/5irr]^)]. 

The time exponent of this final relation is the same 
as Brownian-coagulation mechanism, but the prefactor 
is larger by a factor of 7 than it. This result is quanti- 
tatively consistent with the recent careful experimental 
studies on droplet phase separation in space by Perrot et 
al. p6fl : They demonstrated the i 1 / 3 law with the twice 
larger prefactor than that of Brownian-coagulation mech- 
anism [see Eq. ([!])], using the relation kd = 12$^ M. 
This means 7 1 / 3 ~ 2.6, if we use the relation kd — 6$^, 
which we believe is the correct one ||. In our estima- 
tion 7 1 / 3 = 1.7 - 2.2 at maximum, and the agreement 
is satisfactory at least on a qualitative level. If we take 
into account the deviation of <j) from 4> e |24[] , 7 1 / 3 be- 
comes smaller to be between 1 and 2. We need a more 
analytical approach for the more accurate estimation of 
7- 

Finally, we consider novel hydrodynamic effects unique 
to fluid systems that are not included in the above calcu- 
lation: (i) Direct hydrodynamic interactions: Inter- 
droplet collision and the resulting shape relaxation (or 
the relaxation of curvature) directly induces hydrody- 
namic flow via Eq. (||) and leads directly to another 



3 



collision. This effect likely leads to v ~ a/n and R ~ t 
as first demonstrated by Nikolayev et al. jl6| and also 
by us JO ], This mechanism is likely important only for 
the case of a very high droplet density |TJ|,[l6|]. Since 
this problem is beyond the scope of this paper focusing 
on the growth law of R ~ t 1 / 3 , we do not discuss it in 
detail here, (ii) Coupling between hydrodynamic 
and diffusion modes: The quick hydrodynamic coales- 
cence after a collision could leave the droplet out of local 
equilibrium and temporally cause a strong diffusion field 
around it. In our previous papers [ p7[ , we discuss the 
effects of interface quench induced by the quick hydro- 
dynamic reduction of the interfacial energy, focusing on 
bicontinuous phase separation. We expect a similar, but 
qualitatively different, dynamic effect for the process of 
a quick hydrodynamic coalescence after interdroplet col- 
lision for droplet phase separation 12 l| : The concen- 
tration profile around a droplet just after the collision 
is about the same as that before the collision since the 
diffusion process is retarded from the hydrodynamic one. 
Note that only diffusion can change the concentration 
profile. Further, the radius of the droplet increases from 
R to R' immediately by collision, satisfying the volume 
conservation, (a) This increases the interaction strength 
by a factor of (R' /R) for a droplet experiencing colli- 
sion and increases the probability of the subsequent col- 
lisions. This situation can be realized since the character- 
istic time of droplet fusion tf US i on ~ rjR/a is generally 
much shorter than the characteristic material diffusion 
time for a droplet idi/f ~ R 2 /D^. Using the relation 
a ~ a k B T/£ 2 , we obtain tf usion /t diff ~ £/lQR. The 
strong excess diffusion fields induced by collision lasts 
for about tdiff after the interdroplet collision, (b) More 
importantly, the collision significantly increase the local 
value of V 2 S(j) (or V 2 /x) itself around the droplet for a 
period of tdiff and, thus, produce the strong attractive 
interactions between the droplet experiencin g th e colli- 
sion and the surrounding droplets via Eq. (|l2|). This 
leads to the subsequent collisions. 

All these effects [(f), (h)(a), and (ii) (b)] likely 
strengthen "gradient-induced coupling mechanism" se- 
lectively for droplets experiencing collision and lead 
to the "collision-induced collision" phenomena [ p~2|Jl3| 
that a droplet experiencing collision has more probabil- 
ity of the subsequent collisions. Such behavior is actu- 
ally observed experimentally in droplet phase separation 
[ p2| . This mechanism can be important when t co u < 
tdiff- Here t C oll is the characteristic interval of droplet 
collision by the hydrodynamic translational motion of 
droplets. It is estimated as t coU (L) — f^ 2 dr/v^ ax = 
L 3 1 (24DffR'j), where L is the average interdroplet dis- 
tance. Since L 3 = (4ic/3)R 3 /$ d , t coU = nR? / '187$ d D R . 
The ratio t coU /t d iff is, thus, estimated as t coU /t di ff ~ 
7ri?/187<I> d £ - 0.02(i?/£)$ d . Thus, this effect is likely 
important for the stage where R <~ 200£, for example, 



at $d ~ 0.25. This relation also suggests that in the 
very late stage the evaporation-condensation mechanism 
plays an important role. However, there have so far been 
no reliable experiments on such an extremely late stage. 

In summary, we study the coarsening mechanism of 
droplet phase separation. We find that the interdroplet 
interaction caused by the anisotropic coupling of concen- 
tration fields can accelerate droplet coarsening, in con- 
trast to the conventional understanding that there are 
no interdroplet interactions, which is the precondition 
of Brownian-coagulation mechanism. We obtain the do- 
main coarsening law of R ~ i 1 / 3 , whose prefactor is at 
least larger than that of Brownian-coagulation mecha- 
nism. We call this new mechanism of droplet coarsen- 
ing "gradient-induced coupling mechanism". We also 
demonstrate the possibility that the strong coupling be- 
tween velocity and diffusion fields may lead to a new 
qualitative effect, " collision- induced collision". Further 
theoretical and experimental studies are highly desir- 
able for the clear understanding of this problem, includ- 
ing another effect of "collision-induced collision via flow" 
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